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We theoretically examine the optomechanical interaction between a rotating nanoparticle and an orbital angular 
momentum-carrying optical cavity mode. Specifically, we consider a dielectric nanosphere rotating uniformly in 
a ring-shaped optical potential inside a Fabry-Perot resonator. The motion of the particle is probed by a weak 
angular lattice, created by introducing two additional degenerate Laguerre-Gaussian cavity modes carrying 
equal and opposite orbital angular momenta. We demonstrate that the rotation frequency of the nanoparticle 
is imprinted on the probe optical mode, via the Doppler shift, and thus may be sensed experimentally using 
homodyne detection. We show analytically that the effect of the optical probe on the particle rotation van¬ 
ishes in the regime of linear response, resulting in an accurate frequency measurement. We also numerically 
characterize the degradation of the measurement accuracy when the system is driven in the nonlinear regime. 
Our results are relevant to rotational Doppler velocimetry and to studies of rotational Brownian motion in a 
periodic lattice. 

OCIS codes: (080.4865) Optical vortices; (140.4780) Optical resonators; (260.6042) Singular optics; 

(280.3340) Laser Doppler velocimetry. 
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1. Introduction 

Cavity optomechanics, which realizes the coupling of op¬ 
tical radiation to mechanical motion, is a versatile plat¬ 
form for sensing technologies [Hi- Prime instances of 
this capability are the ultrasensitive linear displacement 
detectors used for gravitational wave interferometry Q 
and atomic force microscopy 0, as well as recently de¬ 
veloped accelerometers Q, magnetometers 0 and ther¬ 
mometers 0. Given its effectiveness in supporting de¬ 
tection functionalities, it is relevant to inquire if cavity 
optomechanics can be used for sensing nanomechanical 
rotation, a capability critical to fields such as nanotech¬ 
nology [1^, Doppler velocimetr\ , atomtronics 

[a and statistical mechanics [isl Il9j| . 

In this article, we consider the rotational motion of a 
nanoparticle inside a Fabry-Perot cavity. While most of 
cavity optomechanics has addressed vibrational motion 
0 , there is some precedence for studies of rotational mo¬ 
tion [2^ . Torsional oscillators have been used to address 
angular deflections [2ll - l^ : but these systems exhibit 
restricted, rather than free, mechanical rotation. Pseu¬ 
dospins such as cold atoms [2^ and Bose-Einstein con¬ 
densates [131 have been coupled to optical cavity fields. 
However, these systems do not actually involve physi¬ 
cal rotation; instead, their internal degrees of freedom 
mathematically map on to rotor models. Optomechan¬ 
ical gyroscopes have been considered but they require 
rotation of the optical cavity itself [23 - l30l| . 

In a pioneering theoretical suggestion, the orbital mo¬ 
tion of a nanoparticle around a spherical whispering 
gallery mode resonator was shown to generate optical 


signatures of rotational motion [3lj |. However, in that 
case the particle is kept in rotation by a torque arising 
from the optical mode itself. This torque thus ultimately 
determines the particle rotation frequency, and its pres¬ 
ence is undesirable for sensing an arbitrary, externally 
imposed nanomechanical rotation rate, or for investigat¬ 
ing dynamics that are driven only by a thermal torque, 
i.e. Brownian motion. Further, the optical torque in the 
earlier proposal is nonconservative, leading to compli¬ 
cated non-Hamiltonian dynamics. Lastly, the nanoparti¬ 
cle in that case can only interact with three values of the 
optical angular momentum, namely 0 and ±/i [^ . This 
constraint places a restriction on the general exploita¬ 
tion of the orbital angular momentum (0AM) degree of 
optical freedom. 

In this article, we investigate a cavity optomechan¬ 
ical scheme which is sensitive to full mechanical rota¬ 
tion and does not require motion of the cavity. Unlike 
previous work, our proposal utilizes two optical modes, 
a trap mode which exerts no torque on the nanoparti¬ 
cle, and a probe mode which exerts a small (conserva¬ 
tive) torque with negligible effect on the particle rotation 
in the regime of linear response of the system. Lastly, 
the particle in our scheme can couple to any value of 
the optical 0AM 1; this opens up the possibility of a 
cavity-based version of standard rotational Doppler ve¬ 
locimetry 0, El, [H [il|. In order to demonstrate our 
scheme, we characterize theoretically below the measure¬ 
ment of the constant rotation rate of a uniformly orbiting 
nanoparticle. 
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2. Configuration 

In this section we describe the physical setup underly¬ 
ing our system of interest. We first briefly address the 
trapping of the nanoparticle in a cavity, and then offer a 
more detailed description of the detection of its rotation. 


2.A. Trap mode 

We consider a Fabry-Perot cavity with two highly re¬ 
flective mirrors, as shown in Fig. [1] Two optical fields 
excite the cavity. The first, trapping, field is an intense 
Laguerre-Gaussian beam LGi' carrying orbital angular 
momentum V 7^ 0, at a wavelength At resonant with the 
cavity. The practicality of exciting such cavity modes 
has been considered earlier, both theoretically [s^ as 
well as experimentally [3^ . For such a mode, the corre¬ 
sponding intracavity intensity is [33| 


Iv{^ 


P' 


|/'|!r7ri?2 



e cos^ {ktz ), (1) 


where r and z are radial and axial cylindrical coordi¬ 
nates, respectively, T is the intensity transmission coef¬ 
ficient of either cavity mirror. 


h = 
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is the wave-vector, is the input power, and 

, 1/2 
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( 2 ) 
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is the mode spot size [^, wo(» \t) being the beam 
waist. 






Fig. 1. (Color online). The optomechanical nanorotation 
sensing scheme proposed in this article. Cavity mirrors Mi 
and M 2 confine an intense Laguerre-Gaussian beam LGi/. 
This annular beam traps a dielectric nanoparticle of mass m 
on a ring of radius R. The dielectric rotates at a constant 
mechanical frequency uJm,s- Its motion is probed by a weak 
angular lattice formed by two additional Laguerre-Gaussian 
beams LG±i. In the figure, I = 2, resulting in four lattice 
sites on the trapping ring. The rotation of the particle is 
detected by homodyning the probe mode LG±i after it exits 
the cavity. 


The energy of interaction between the trapping optical 
mode and a dielectric sphere of mass m and radius ra < 
At (Rayleigh regime) in the dipole approximation is 

20 - 

= - h'if), (4) 

ceo 

where 

“ = (^) • <"> 

is the dielectric polarizability, eg the permittivity of vac¬ 
uum, er the relative dielectric constant, c the velocity of 
light and 


V = 


47rrg 

3 


( 6 ) 


the volume of the dielectric sphere. Here we have as¬ 
sumed that the particle is small enough that it does not 
break the rotational symmetry about the axis of the cav¬ 
ity sufficiently to resolve the several degenerate Hermite- 
Gaussian modes that superpose to yield the trapping 
Laguerre-Gaussian mode. 

The interaction energy is minimized at r = i?, z = 0; 
the dielectric particle is therefore trapped on this ring, 
which corresponds to a field maximum. The presence of 
the particle at the antinode of the optical trapping field 
will shift the cavity resonance. However, this shift can 
be easily compensated by changing the trap laser fre¬ 
quency by a small amount. A large P'^ is chosen such 
that the trapping along the radial (r) and axial (z) direc¬ 
tions is stiff. The corresponding particle displacements 
are then small and effectively decouple from the par¬ 
ticle’s azimuthal motion Furthermore, such small 
displacements may only pull the trap mode frequency 
by an amount which is much smaller than the cavity 
linewidth, and is thus negligible Q. In the remainder 
of the article, we therefore ignore the axial and radial 
motion of the dielectric. Further, the intense trapping 
optical field will be considered to be a parametric quan¬ 
tity and not a dynamical variable. 

The azimuthal motion of the particle is free, and stable 
if the rotational energy of the particle does not exceed 
the depth of the confining radial well. We will assume 
that the dielectric particle experiences a constant torque 
T causing it to rotate around the cavity axis. This torque 
can be enforced by magnetic actuators if the particle is 
chosen to be a paramagnetic bead [s^, or if its surface 
is partially coated with metal, making it a “dot Janus” 
particle |36j | . Such rotation techniques have been demon¬ 
strated experimentally, albeit in liquid media; however, 
several experiments have also recently trapped and ro¬ 
tated microparticles in free space [H, • Rather than 

consider in further detail the mechanism of rotation, 
which may be implemented in a number of ways, we 
proceed to focus instead on our detection scheme. 
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2.B. Probe mode 

The rotation of the particle is weakly probed by a sec¬ 
ond optical field of wavelength Ap 7 ^ At, near-resonant 
with the cavity, and consisting of a superposition of two 
degenerate counter-rotating Laguerre-Gaussian beams 
LG±i with orbital angular momenta -1-1 and —I respec¬ 
tively. Near the trap ping ring, the corresponding mode 
function obeys 

. 

where loq ^ Ap = 27r/fcp. In practice, the polarizations 
of the trapping (At) and probe (Ap) fields may be chosen 
orthogonal, in order to eliminate interference effects. As 
can be seen from Eq. the effect of the LG±i beams 
is to create an angular lattice with 21 sites on the trap¬ 
ping ring (Fig. [T] displays the case I = 2). This lattice 
breaks the continuous azimuthal symmetry about the 
cavity axis and can thus sense the rotation of the parti¬ 
cle [IJ, [I 3 , as we show below. 


the rotor, we use the classical angular momentum Lz and 
the explicitly periodic angular displacement variable 

Ui = ( 11 ) 


where (j) takes on continuous values. Using the variables 
just described, the coupled-mode equations of motion in 
the frame of the laser driving the cavity with the LG± 
fields are 


I* A'-^([/z + t/n -2|a + V7azn,(12) 


Ui = 


i2lUiLz 


(13) 


Lz = -ImLz - 2ilhg(l) {Ui - U*) \a{t)f -I-r -|- 


We now describe each of the Eas. (fT^ - (IT4l) . In Eq. (IT^ . 
which has been derived using the formalism of Haus [39j , 



(15) 


2.C. Optorotational coupling 

We now characterize the light-matter coupling for our 
system, using a classical model. A fully quantum me¬ 
chanical treatment will be the subject of later work. We 
denote the cavity resonance frequency of the probe LG±i 
modes in the absence of the dielectric as ujc- In the pres¬ 
ence of the dielectric, this frequency shifts in a way that 
depends on the angular position (j) of the dielectric [ 2 ^ . 
Since < {Xp,L), where L is the length of the cavity, 
the shifted frequency uj{(p) can be evaluated using the 
perturbation-theoretic Bethe-Schwinger formula (38| 

^ , _ /v(er - 1) |V>;(r)|^dr 

Wc dr 

where 14 is the cavity mode volume. Using Eq. 0 in 
the formula m gives 

uj{(j)) = ujc — g{l) cos^ I (j), (9) 


where A = ujd — ujc is the detuning between the fre¬ 
quency LOd of the LG±i driving laser and Wc, 7 is the 
cavity loss rate, and am = \/PmlhuJc., where Pm is the 
input power. Classical laser noise has been neglected, 
as well as the thermal noise contribution to the radia¬ 
tion mode, which is very small at optical frequencies. 
Vacuum fluctuations in the optical mode have also been 
ignored (but see below). Comparing Eq. (IT^ to the 
analogous equation for vibrational optomechanics im¬ 
mediately clarifies that the phase of the optical field is 
sensitive to the (Ui P U*)/2 = cos2l(j) quadrature of the 
mechanical rotation, and that this sensitivity improves 
linearly with the cavity finesse Q . 

In Eq. m, I = mP? is the dielectric particle’s mo¬ 
ment of inertia about the axis of rotation, and we have 
implicitly used the relation (j) = Lzjl when taking the 
time derivative of C/;. 

Equation m describes a rotor damped at a rate 7 ^., 
and Tin is a Langevin torque with zero mean and the 
two-time fluctuation correlation 


with 


9{l) = (er 


1 ) 





( 10 ) 

where r[(l-|-1)/2] is a Gamma function. We note that in 
contrast to previous work the dielectric particle is 
far away from the cavity mirrors, and thus no apprecia¬ 
ble rearrangement of charges occurs on the cavity walls, 
and the standard form of the optomechanical coupling 
applies. 


3. Equations of motion 

To describe the dynamics of the optical field, we use the 
classical variable a{t), such that |a(t)|^ is the number of 
photons in the cavity at time t. In order to characterize 


{5Tin{t)5nnit')) = 2IjrnkBT6(t - t'), (16) 

signifying white noise. Unless mentioned otherwise, 
we will assume that the externally applied torque is 
larger than the torque due to the optical lattice, i.e. 
T>2ilhg{l){Ui-U^)\a{t f . 


3.A. Steady state 

The steady-state of the system can be found readily by 
equating the time derivatives in Eqs. (HDl-dll to zero, 


_ Lz,s _ T 

I 


(17) 


Uls = 0 , 


(18) 
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^/lain 


+ ( 2 ) 



(19) 


Since the particle is in a non-equilibrium steady state, 
Equations (dZli-dlll) are to be interpreted as statements 
about time averages. For example Eq. (HI can be re¬ 
written as 


Ui,s = m, 

= (cos 2l(j) + i sin 2Z(/)) 

= {cos2l(j)) + i{sm2l(j)), (20) 


where the brackets denote an average over one period of 
mechanical rotation. Thus, while can never equal 
zero as a function, it can vanish on average, since the 
cosine and sine average to zero during one rotation. 

Importantly, from Equation (II3 it can be seen that 
the torque due to the probe lattice makes no contribu¬ 
tion to the average rotation rate w^.s- This is because 
the lattice accelerates the particle over half of its mo¬ 
tion while decelerating it (almost) equally over the other 
half. This suggests that a rotation sensing (velocimetry) 
scheme based on our proposal can circumvent, at least 
to lowest order, the effect of the optical probe on the 
mechanical motion, and thus provide an accurate mea¬ 
surement of the undisturbed rotation frequency. This 
statement will be supported both analytically as well as 
numerically below. 

It is also worth noting that working directly with (j) = 
Lz/1 instead of Eq. (fT51) would have implied no rotation 
in the steady state, which is unphysical, since we have 
assumed that the external torque r always overcomes 
the weak angular lattice. 

3.B. Linear response 

We now examine the linear response of the system, by 
writing each dynamical variable as a sum of its steady 
state value and a fluctuation Q, e.g. 

a = Os -|- 5a, (21) 

Using Eqs. (fT^ - dT^ and retaining only terms linear in 
the fluctuations, we find 

5a = -i^as {5Ui + 5U!) + (f A' - 5a, 

72/ 

5Ui = —Lz,s5Ui, ( 22 ) 

5Lz = -^m5Lz + 2ilhg{l) |as|^ {5Ui - 5U*) + 


This set of linear first-order differential equations can 
be readily solved analytically for any initial conditions 
5a{0),5Ui{0) and 5Lz{0). The full solution is rather in¬ 
volved; here we consider only a revealing limit: on res¬ 
onance (A' = 0), and ignoring cavity losses (7 = 0), we 


find the Fourier transform of 5a(t) is 

5a {uj) = A5 (oj) + B*5 (oj — cvg) — B5 (w -I- uJs) , (23) 

where 

A = V^5a{0) + B - B*, (24) 

^ ^ ^/^asg{l)5Ul{0) ^^5) 

^m,s 

and 

UJs — ‘2luJm,s- (^^) 

The Dirac delta at w = 0 on the right hand side of 
Eq. (1^ corresponds to the optical frequency, since we 
are in the frame rotating at the cavity resonance. Im¬ 
portantly, the Dirac deltas at uj = ±Ws correspond to 
sidebands which are created as the particle rotating at 
the frequency uJm,s encounters 21 optical lattice sites. 
These sidebands fundamentally arise from the rotational 
Doppler shift [13,113 imprinted on the cavity photons by 
the mechanical motion, and are analogous to the side¬ 
bands at the harmonic oscillator frequency seen in stan¬ 
dard vibrational optomechanics experiments Q. The 
rotation frequency peak at cug can thus be recovered by 
homodyning the probe beam LG±i once it has exited 
the cavity [3|. Interestingly, I maybe treated as an ad¬ 
justable parameter by the experimentalist to allow con¬ 
venient detection for a given rate of rotation. 

Equation (1^5)1 shows that the linear response is sensi¬ 
tive to the average mechanical rotation rate ujm,s, and is 
insensitive to the action of the probe field on the mechan¬ 
ical motion. It can readily be verified that this conclu¬ 
sion holds true even in the presence of detuning (A' 0) 

and optical damping (7 7^ 0). 

3.C. Dynamics 

To analyze more comprehensively the dynamics of the 
system, we now numerically simulate the full nonlinear 
model [Eqs. (fT^ - (fTi)) ]. 

3.C.I. Simulation parameters 

To make the simulation realistic, we include vacuum op- 


tical noise by using 



V7am y/jain + y/^5ain, 

(27) 

in Eq. 

([T^. with Q 



{5ain{t)) = 0, 

(28) 

and 

{5ain{t)5a\^{t')) = 5{t-t'), 

(29) 


All stochastic processes are modeled as Gaussian white 
noises, and each simulation run includes a single realiza¬ 
tion of the noise. 
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For the cavity, we assume L = 25 cm, and T = 
8 X 10“"^, which lead to 7 ~ 150 KHz. The radius of 
curvature of both mirrors is assumed to be 40m. For 
the trapping field, At = 1064 nm, leading to ojq = 2 
mm, and V — 2, implying i? ~ 2 mm; also, we as¬ 
sume = 180 mW. For the probe field, we assume 
A = 980 nm. A' = 0,1 = 2, yielding ujc — 10^® Hz. 
The dielectric particle is chosen to be a polystyrene 
(cr = 2.5) sphere of radius = 150 nm. This gives 
a = 1.2 X 10“^^m^ and g{l = 2) ~ 25 mHz. Us¬ 
ing the polystyrene density p = 1050 kg m“^, we find 
m = 1.5 X 10“^® kg, and I = 62.5 x 10“^^ kg m®. From 
earlier work, 7p ~ 60 Hz at a background pressure of 
^ 10 mBar [^; also, we assume r = 2.5 fN/rm, which 
is within demonstrated capabilities [s^. These assump¬ 
tions lead to Lz^s = 8.3 x 10“^^kg m®s“^ ~ 10^°fi, justi¬ 
fying our classical treatment of the rotor, and implying 
a rotational kinetic energy of L® g/2I ~ 107 K. This is 
safely lower than the radial trapping well depth ~ 360 
K; thus the particle remains bound as it rotates. The 
optical intensity experienced by the particle is about 120 
kW cm“®, several orders of magnitude below the thresh¬ 
old for heating due to absorption [d^. Since <C Xp, 
the effect of recoil heating on the particle motion is neg¬ 
ligible [13,113 ; also we find that the scattering-induced 
photon loss rate from the cavity is small 30Hz) [d^ 
compared to 7 calculated above and does not greatly 
affect the cavity finesse. 

3.C.2. Dynamics: linear regime 

Using the standard methods of dynamical systems the¬ 
ory, we have ensured that all the configurations discussed 
in this article are mechanically stable [dj. Starting 
from a cavity empty of probe radiation, and a parti¬ 
cle at rest at (p = 37°, we now consider the behavior 
of the system for various probe powers Pin- In Fig. [5] 
(a)-(d), Pin = 2/rW, which lies in the regime of lin¬ 
ear response, but is adequate for carrying out homo¬ 
dyne detection. Figure (a) shows the number of probe 
photons in the cavity reaching an equilibrium value; (b) 
shows that the imaginary part of the cavity field, which 
can be measured using homodyne detection, displays si¬ 
nusoidal oscillations due to mechanical modulation; (c) 
shows the Fourier transform of (b) with a frequency peak 
at uifi = 72.463761 ± 10“®Hz, which accurately reflects 
the nominal value ujs from Eg. (1261) (d) shows the time 
dependent mechanical rotation rate uim(t) rising to the 
steady state value uJm,s = w^/d. Further numerical ex¬ 
ploration shows that the limits on the frequency reso¬ 
lution arise mostly from thermal and to a lesser extent 
from photonic noise. 

3.C.3. Dynamics: nonlinear regime 

In Figs.[5](e)-(h), Pin = 10 mW, which lies in the regime 
of the system’s nonlinear response, i.e. when the vari¬ 
ables deviate appreciably from their steady state values. 
Figure (e) is similar to (a); (f) however shows a more 
complicated time evolution than (b) and suggests the 
presence of multiple underlying Fourier frequencies; in¬ 





Fig. 2. Dynamics of the proposed optomechanical system for 
T = 25mK (other parameters are stated in the text) and 
various regimes of the probe power: (a)-(d) Pin = 2 /rW, (e)- 
(h) Pin = 10 mW. (i) Homodyne peak location versus probe 
power, (j) Potential energy of the optical lattice. The plots 
are described in detail in the text. 


deed, (g) shows a peak at = 55.26934 ± 0.00002Hz, 
and higher sidebands at integer multiples of a;)^; (h) 
shows large optical modulation of the mechanical rota¬ 
tion rate and may be compared to (d). 

3.C.4. Velocimetry 

The two system response regimes, linear as well as non¬ 
linear, are summarized together in Fig. Hi), which shows 
the homodyne-detected frequency uJd versus the probe 
mode power Pin- In the linear response regime of low 
probe power, as Pin is varied, the homodyne signal fre¬ 
quency UJd stays constant and close to the nominal value 
ujs - The identification of this regime of accurate response 
is the central result of this article. In this regime, the 
potential energy of the lattice [Fig. [2Dj)] is much smaller 
than the kinetic energy of the particle, and modifies the 
rotation frequency very little from the nominal value 
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^m,s- While Figs. (I2])(a)-(d) display the linear regime 
results for T = 25 mK, runs for T = 300 K [not shown 
in Fig. ([2])] yield LOh = 72.4738 ± 10“^ Hz. These mea¬ 
surements are of at least an order of magnitude higher 
resolution than state-of-the-art velocimetry experiments 
[ilElIlliaiii, implying that cavity-based rotation 
sensing may be a direction worth exploring for rotational 
Doppler velocimetry. While we have investigated the 
detection of low angular velocities (few tens of Hz) for 
simplicity, in practice it might be more convenient to 
work with faster rotation rates, due to the presence of 
low frequency noise as well as fluctuations in background 
pressure and the applied torque r. 

In the nonlinear regime of high probe power [beyond 
the critical power , supplied as a guide to the eye, in 
Fig. ©(i)], the particle kinetic and lattice potential en¬ 
ergy are comparable, and the negotiation of the lattice 
adds a large delay to the orbital period of the parti¬ 
cle. This lowers the average particle rotation rate sig¬ 
nificantly from its unperturbed value uJm,s, which can 
no longer be measured accurately. Eventually, for very 
large Pin, the particle is trapped in one of the lattice 
minima. 

4. Conclusion 

We have theoretically explored the optomechanical cou¬ 
pling of nanorotation to a cavity mode carrying orbital 
angular momentum. Our results suggest that optical 
cavities deserve to be investigated further as a means of 
advancing the technique of rotational Doppler velocime¬ 
try. 

Our analysis can be extended to the case r = 0, in 
which case the only external torque present is thermal, 
and the dielectric particle travels on the angular peri¬ 
odic lattice created by the probe mode. This would 
correspond to Brownian motion in a periodic angular 
lattice, which is of interest, for example, to studies of 
dielectric relaxation This would likely require very 
low background pressures. It would also be interest¬ 
ing to generalize our classical treatment to the quantum 
regime, to determine, for example, what the minimum 
detectable rotation rate - i.e. mechanical angular mo¬ 
mentum - is. Lastly, our analysis may be extended to 
describe the motion of an atom in place of the dielectric, 
as the optomechanical coupling will be of the same func¬ 
tional form in the angular coordinate, if the probe laser 
is detuned from an atomic resonance, and the interaction 
is consequently due to a dipole potential. 

We are grateful to S. Agarwal, B. Zwickl, M. Ven- 
galattore, H. Shi, B. Ek and B. Rodenburg for useful 
discussions, and to the Research Corporation for Science 
Advancement for support. 
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